Spectral determinant on quantum graphs 

Eric Akkermans^'*'"'" , Alain Comtet*, Jean Desbois*, 
Gilles Montambaux* and Christophe Texier^* 

1st February 2008 



t Department of Physics. Technion, Israel Institute of Technology, 32000 Haifa, Israel. 

* Laboratoire de Physique des Solides. Universite Paris-Sud, Bat. 510, F-91405 Orsay Cedex, France. 

* Laboratoire de Physique Thcorique et Modeles Statistiques. Universite Paris-Sud, Bat. 100, F-91405 
Orsay Cedex, France. 

* Departement de Physique Theorique, Universite de Geneve, 24 quai Ernest Ansermet, CH-1211 
Geneve 4, Switzerland. 

Abstract 

We study the spectral determinant of the Laplacian on finite graphs characterized by their 
number of vertices V and of bonds B. We present a path integral derivation which leads to 
two equivalent expressions of the spectral determinant of the Laplacian either in terms of a 
V x V vertex matrix or a 2B x 2B link matrix that couples the arcs (oriented bonds) together. 
This latter expression allows us to rewrite the spectral determinant as an infinite product of 
contributions of periodic orbits on the graph. We also present a diagrammatic method that 
permits us to write the spectral determinant in terms of a finite number of periodic orbit 
contributions. These results are generalized to the case of graphs in a magnetic field. Several 
examples illustrating this formalism are presented and its application to the thermodynamic 
and transport properties of weakly disordered and coherent mesoscopic networks is discussed. 

1 Introduction and main results 

This work is devoted to the study of the spectral properties of the Laplacian operator on finite 
graphs. This problem has already a long history. The properties of the Laplacian describing 
free electrons on networks made of one-dimensional wires have been investigated in the context 
of organic molecules |IJ . Subsequently, this approach has proved useful to study superconduct- 
ing networks using linearized Ginzburg-Landau equations [^J, the vibration properties of fractal 
structures such as the Sierpinski gasket ||, the adiabatic conductances of a network in an in- 
homogeneous magnetic field [Q, ||] or the behaviour of disordered systems in a magnetic field 



|, y, 1C, 11]. More recently, it has been shown JUj, 13, 14] that the Laplacian on graphs pro- 
vides an interesting framework to study the onset of quantum manifestations of chaos. In these 
examples, the quantity of main interest is the spectrum of eigenenergies of the Schrodinger equa- 
tion, or the diffusion equation, defined for each wire of the network with appropriate boundary 
conditions at the vertices. Finally we should also mention the relevance of the spectral determi- 
nant of the Laplacian on graphs in the context of statistical field theory and more precisely for 
the problem of the triangulation of random surfaces. There, the corresponding homology and 
(dual) cohomology groups provide the corresponding graphs over which the Laplacian is defined 

The problem of the spectrum of graphs has been also investigated thoroughly in the mathe- 
matical literature (l|, |l7|, H, ||, |0|. The more specific question of deriving a trace formula for 
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the heat kernel (partition function) of the Laplacian operator on a graph has been studied by 
Roth |||, |||. 



Let us start by recalling some already known results. The partition function Z(t) of the 
scalar Laplacian operator A with appropriate boundary conditions, is defined as 

Z(t) = Tv{e tA } . (1) 

Evaluating the trace over a set of eigenstates of —A of eigenvalues E n , we obtain 

Z(t)=Xy***- (2) 

n 

The spectral determinant S(j) is formally defined by S{^) = det(— A + 7) = Iln(7 + E n ) up to 
a regularization independent of 7. It is such that 



poo a 

/ dtZ(t)e-^ = — In 5( 7 ) 
Jo 07 



(3) 



On a graph in the presence of a magnetic field, a compact form for the determinant of the 
operator 7 — (d x — iA(x)) 2 has recently been obtained 11]: 



S{l) =7^ II My/lLp) det(M) , (4) 

(a/3) 

where V is the number of vertices, B is the number of bonds and l a p is the length of the bond 
(a/3). M is a V x V matrix whose elements are [^, ||, [/], |i~lfl : 

M aa = ^coth(V7Wi) ( 5 ) 



i=l 



M a « = — — — - if (a/3) is a bond (6) 

= otherwise , 

where the summation is taken over the m a neighbouring sites of the vertex a. 9p a = J A(x)dx 
is a phase equal to the circulation of the vector potential A(x) along the bonds (af3). This ex- 
pression has been obtained by a calculation of the Green's function on the graph. The derivation 
is recalled in appendix [A|. 

Although it is well known that the spectrum of the Laplacian on a graph is given by the zeros 
of det(M (7 = -E)) = 0, |, §], formula @ is richer since it contains non trivial multiplicative 
factors which depend on 7. If the relation of 5(7) to the matrix M(j) allows a very simple 
determination of the spectrum, it cannot help to understand the structure of the spectrum in 
terms of periodic orbits on the graph. To do so, it is more convenient to introduce a "link 
matrix" E whose elements Er^i^s relate oriented bonds (a/3) and (ixv). This 2B x 2B matrix 
is defined as[]: 

= ^e-v^W+i^ if/? = /U (7 ) 

JL _ l^j e ~V7 W+ i0 c</3 if / P = V (8) 
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a = v 



otherwise . (9) 



1 The matrix E is the same as the matrix Sb in 
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The spectrum is now given by det(l — E(^y)) = 0. The trace expansion of det(l — E) leads to a 
description of the spectrum in terms of periodic orbits |R| . 

Finally, a trace formula has been obtained by Roth, which expresses the partition function 
as a contribution of such periodic orbits: 

L V - B 1 „ - i(c) 2 

ZW = 57S + — + vsS™>e— (io) 

where the sum runs over closed orbits C on the graph. C is the primitive orbit associated with 
C . The geometrical factor a(C) will be defined in section |6| 

We now give an overview of the paper and describe rapidly our main results. In section |2| 
we recall briefly how transport and thermodynamic properties of weakly disordered and coher- 
ent conductors can be related to the spectral determinant || ^, pi]]. In particular the weak- 
localization correction to the conductance is a physical measure of the spectral determinant. 

After having set up the notations and definitions in section ||, we present in section || a 
path-integral formulation for the calculation of the spectral determinant of the Laplacian on a 
graph. It allows us to obtain an expression of the spectral determinant in terms of the vertex 
matrix (||]6|). The interest of this approach relies on the fact that the path-integral on each 
bond involves the propagator of a two-dimensional harmonic oscillator for which the role of the 
time is played by the length along the bond. This problem corresponds to a zero-dimensional 
Gaussian field theory, so that the path-integral is easily calculated using standard quantum 
mechanics. This approach has many advantages like to allow a simple generalization to other 
types of boundary conditions, or to permit very easily the elimination of vertices of coordinence 
2 (appendix O). 

In section |5|, we derive a dual and equivalent expression of the spectral determinant in terms 
of the link matrix (0,||): 

S{ 1 )= 1 X ^e^ L (j{m a )2- B det(l-£( 7 )). (11) 

a 

By a trace expansion of the determinant det(l — E), we find in section || the following result 
S( 7 ) = 1 Y ^e^ L (j\m a )2- B \{ (l - a(C)e-^ l ^ +w ^) , (12) 

a q 

which allows us to recover the Roth trace formula for the the partition function ( |l0| ) , demonstrat- 
ing the equivalence between (||) and (|io|). The trace formula (^) for the spectral determinant 
involves the contribution of an infinite number of periodic orbits C and bears some similarity 
with the Selberg trace formula |2^] or the semi-classical trace formulae for chaotic Hamiltonian 
systems ]24|]. In section ^ we show that the spectral determinant involves the contributions of 
a finite number of periodic orbits only. This provides a diagrammatic method to compute the 
determinant. This method is applied to a particular example. 

In section ^ the path integral formulation is shown to be a good starting point to generalize 
these results to the case where a magnetic field is applied to the network. In this way we obtain 
a generalization of the Roth formula ( |i~0| ) . Sections ^ and |l(] discuss how to extend the previous 
results to the case of more general boundary conditions, and provide a discussion of the low 
energy behaviour of the spectrum and the existence of zero modes, respectively. 



Section |11| discusses the scattering problem when the graph is connected to an infinite lead. 
In this case the relevant information is contained in a phase shift which is shown to be related 
to the ratio of two spectral determinants calculated with different boundary conditions. Finally, 
section ^ gives an illustration of the formalism for the case of the complete graph K n . 
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2 Phase coherent properties of disordered conductors 



In this section, we briefly recall how transport and thermodynamic properties of weakly disor- 
dered and coherent conductors can be related to the partition function Z(t) or equivalently to 
the spectral determinant £(7). A more detailed derivation can be found in references |8|, [| |l0| 
and references therein. These expressions are quite general and are not specific to graphs. But 
their calculation on any network made of quasi-one-dimensional diffusive wires is straightforward 
with the expression (B) [Tlfl . 

The first step is to write the physical quantities in terms of the classical return probability 
P(r, f, t) solution of the diffusion equation (we set h = c = 1) 



d / - -a 2 



P(r,f',t) = 5{f-r') (13) 



where — e is the electron charge, D is the diffusion coefficient. In the other sections of this 
paper, we have set — 2e = 1 and D = 1. The return probability has actually two contributions, 
a purely classical P c i which is solution of Eq. (^) with A = and another one Pj n j which is the 
result of constructive interferences between electronic propagators associated to time-reversed 
trajectories. This second contribution only exists when the system is phase coherent. In the 
presence of a magnetic field, time reversed trajectories accumulate opposite phases. This is why 



Pint obeys Eq. ( |13|) with an effective charge 2e |2jJ. In this section, we use a unique notation 
for the two contributions. 

The solution of the diffusion equation ( |l3|) has the form 

P(?y,t) =e u (t)^2M^n(r')e- Ent (14) 

n 

where 9n(t) is the Heaviside function. The eigenvalues E n and the eigenfunctions ip n are solutions 
of 

- D (V + 2iel) 2 Mr) = ^n^n(r) • (15) 

The integral over the whole space of the return probability is precisely the partition function of 
the Laplace operator —A, or —(V + 2ieA) 2 in a magnetic field, namely 

Z(t) = I drP(r, f,t) =Y^ e ~ Ent ■ ( 16 ) 

n 

As recalled in the introduction (Eq. (||)), the Laplace transform of Z(t) is related to the spectral 
determinant S(-y). 

The electrical conductivity is a current-current correlation function and expresses in terms of 
products of two electron propagators. These propagators have uncorrelated phase factors which 
cancel after disorder averaging. This corresponds to the classical Drude-Boltzmann conductivity 
o"o = e 2 Dpo, where po is the density of states at Fermi energy. However, pairs of time reversed 
trajectories have the same action and thus the same phase. This constructive interference leads 



to an additional contribution to the conductivity, called the weak- localization correction [25 
which is proportional to the probability of having pairs of time reversed trajectories, i.e. to the 
probability for a diffusive particle to return to the origin. The correction can be written as 

9<= 2 D r°° 

(Aa) = (a)-a = -—- dt Z(t)e^ (17) 

7T I I n 
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with f2 being the volume of the system. The exponential damping at large time is due to the 
breakdown of phase coherence because of inelastic processes. This breakdown is phenomeno- 
logically described by a characteristic time = I/7, related to the phase coherence length^] 
= ^/Dt^,. A magnetic field or a Aharonov-Bohm flux, by breaking the time-reversal symme- 
try, destroys the weak-localization correction. 

In a mesoscopic system, i.e. with a typical size L smaller than L^, the conductivity is known 
to exhibit large variations from sample to sample, with a variance which is universal in the limit 
L/Ls — > 0. The structure of this variance results from different phase correlations between four 
electron propagators (two for each conductance). It can be shown that this variance can also be 
written in the following compact form Q: 

<^ 2 > = ^7vr/ d«Z(t)e-* (18) 



/3vr 2 n 2 

where /3 = 1 if there is time reversal symmetry and (5 = 2 in the absence of such symmetry. 

The magnetization due to the orbital motion of the non-interacting electron gas is known as 
the Landau magnetization. Taking into account electron-electron interactions in the Hartree- 
Fock picture gives an additional contribution, known as the Aslamasov-Larkin correction p6| . 
The structure of this correction implies the product of four wave functions, i.e. two local density 
of states whose disorder averaged product can be related to the Fourier transform of the return 
probability P(f,f,t) ^7j. After spatial integration, one gets: 

where Ao is the dimensionless interaction parameter^. B = V x A is the magnetic field. 

Finally the typical magnetization M typ , defined as M 2 yp = (M 2 ) — (M) 2 , can also be written 
in terms of field derivatives of Z(t). Neglecting electron-electron interactions, the magnetization 
is proportional to the density of states. Thus the typical magnetization involves the two-point 



correlation function of the density of states which has been shown by Argaman et al [3C] to be 
related to the Fourier transform of tZ(t). One obtains: 

M L = ^l dt[Z"(t,B)-Z"(t,0)}— (20) 

where Z"(t,B) = d 2 Z(t,B)/dB 2 . 

Using standard properties of Laplace transforms, the above time integrals can be written as 
integrals of the logarithm of the spectral determinant, so that the physical quantities described 
above can be written as |Q ]: 

i2,= 4 n 2 f) 2 

(^) = - WW < 22 > 

1 f 00 d 2 
K P = p I d 7 i(7- 7i)^ lnS(<n)|J (23) 

(M ee ) = -° / d 7 i^mS( 7 i). (24) 



7T 
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2 Integrals ( [L7|jlqJl9| , po[ ) may not converge at small time. In this case the lower bound of the integral should 
be the elastic time 77, , the smallest time for diffusion. 

3 Considering higher corrections in the Cooper channel leads to a ladder summation |2^, p8| , so that Ao should 
be replaced by X(t) — Ao/(l + Aoln(eFi)) = l/ln(Tot) where To is defined as To = €fG ■ The authors of 
reference H found X(i) = l/ln(e F t). 
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When the integrals do not converge at the upper limit, this limit should be taken as l/r e where 
r e is the elastic time. 



3 Definitions and notation 

Consider a graph Q which consists of a set of V vertices linked by B bonds. Its adjacency 
matrix a a p is a square matrix of size V where a a /3 = 1 if a bond joins the vertices a and (5 
and otherwise (in particular a aa = 0)0. The coordinence of vertex a is m a = X^/3 a a/9- Each 
bond (a/3) of length l a g is identified with a finite interval [0, l a p\ £ R. The total length of Q 
is L = I] bonds ( Q( g) = I Eo,/3 a a/3^/3- We denote by x a/ g the coordinate on the bond (a/3), 
starting from vertex a (it follows that xp a = l a p — x a p). A scalar function on Q is a set of B 
components functions ip( a p){x a p). When there is no possible ambiguity we will neglect to label 
components. The Laplacian on Q acts as the usual one-dimensional Laplace operator on each 
bond: 

d 2 

(^)(a/3) = -r-2-4>(a/3){Xaf3) ■ (25) 
ax a/3 

Its domain is the set of functions that satisfy the following conditions at the vertices 

(i) continuity 

VVft) = 0) = for i = l,---,m a (26) 

(ii) a second condition sufficient to insure current conservation 



E 



d x a ^(al3 t ){x a p z = 0) = (27) 



where d x = and the /?j's label the m a neighbouring vertices of vertex a. The scalar product 
is defined as 

flap 

Wm = Zj te^tafiWWafiW- (28) 

(aP) J ° 

In the search of the eigenvalues of the Laplacian, one may construct the wave function of 
energy E = k 2 in a way that conditions (i) is satisfied: V(a/3) = s inM p (V'a sm fc(^a/3 — ^0/3) + 
V'/? sin/cxo,/?). Imposing conditions (ii) lead to the system of V equations M a p{^ = —k 2 )ipp = 
0, where M is given by (^^|). The eigenvalues of the Laplacian are solutions of det(M(7 = 

-£))=o|, 1, in. 



It is useful to introduce some additional notions. An arc (a/3) is defined as the oriented bond 
from a to (3. Each bond is therefore associated with 2 arcs. A path on Q is an ordered sequence 
of arcs such that the end of each arc coincides with the beginning of the one that immediately 
follows in the sequence. Closed paths will play a special role. The set of closed paths that 
only differ by a cyclic permutation of the arcs will be called a circuit (or an orbit). Among all 
possible circuits we will distinguish those that are primitive. A circuit is said to be primitive if 
it cannot be decomposed as a repetition of any smaller circuit. The length of a circuit C will be 
designated by 1(C). 

As an example, let us consider the graph of figure |l|. This graph contains 6 arcs: 1,2,3 and 
the reversed arcs 1, 2, 3. (1,1,2,2) and (2,2,1,1) are two different paths, two possible repre- 
sentations of the same circuit. (1, 1,1,1,1,1) is not a primitive circuit but (1, 1, 1,1, 1,1,2,2) is 
indeed primitive. 

4 If a given pair of vertices is linked by n > 1 bonds, it is always possible to introduce additional vertices 
(see figure without changing the nature of the graph, to go back to a situation where there is only one bond 
between two vertices. However, the formalism that is presented in this paper can be easily generalized to avoid 
this procedure and minimize the number of vertices by suppressing all vertices with coordinence 2. It is the 
purpose of appendix |c] to discuss this point. 
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Figure 1: A graph with 4 vertices, 3 bonds and 6 arcs. 



4 A path integral derivation of the spectral determinant 

An expression of the spectral determinant S("f) = det(— A + 7) in terms of the determinant of 
a finite matrix has been derived in [11, |l0| (c/. Appendix |A|). The purpose of this section is to 
provide a more direct derivation of this result using a path integral formulation. The spectral 
determinant may be written as 



5( 7 ) = det(- A + 7) = ( ! Vcf>V^ e"^ W> ^~ A +^A 

\J <b on Graph / 



(29) 



> on Graph 

where is a complex field defined on the graph Q. The path integral is performed over all fields 
satisfying the conditions (26,27) and the integral in the exponential along all branches of Q. 
Since the field is continuous at each vertex the path integral may be decomposed on each bond 
{oi(5) as an integration over fields that take fixed values (f> a and <fip on the two sides of the bond 



TT dcj) a d4> a TT / p^e-sJo d ^)(- d *+7)^) (3 ) 

.-.SIZL L : i •/</>(())=<*„ 



vertices a bonds 
(a/3) 



where d(/)d<j> = dRe^dlm^. This involves, after an integration by parts, the following quantity^ 

(a/3) ^*(°)=^ 

= «p(~EfcE^^aA = °)) II r (g ^We4^^l^l^l 2 )(31) 



a=l i=l (a/3) 



For the sum of boundary terms in the exponential we have replaced the sum over the bonds by 
a sum over the vertices, a trick that will be frequently used throughout the rest of this paper. 
Equation (|27| ) implies that the boundary terms vanish. The above path integral, involving a 
zero-dimensional Gaussian field theory, can easily be integrated out using standard quantum me- 
chanical techniques. On each bond the path integral can be expressed in terms of the propagator 
of a two-dimensional harmonic oscillator of frequency uj = «/7 

K 1 Pwe4/»^WIW) = ctr^ifaja) = <&|e-V(-^)|&> , 

J<j ) {Q)=4 >a °" 

(32) 

where 4> = (Re (j>, Im cj>) and M r is a constant independent of 7 depending on the precise nor- 
malization chosen in the definition of the path integral, or in other terms on the choice of 

5 The notation X^(a^) ' ' ' ano1 11(0,3) ' ' ' means, throughout this paper, a summation or a product over bonds 
(and not arcs). 
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regularization made to define the determinant. In the following we will drop this inessential 
normalization constant. Using the expression of the propagator 

G£(£ $') = JU-nkW»)(W a )-««') (33) 

2ir sh(u;x) 

one may re-scale the fields <p a to extract a 7-dependent factor: 

r v 

^(t)- 1 = 7^ / n n G^ lag ($ p j a ) , (34) 

J a=l ( a /3) 

where G x (<f>,4>') = G% =1 ((f), <j) ') ■ This expression will be used as a starting point for the derivation 
of the following section. Using the V x V matrix M introduced in section |l] 

M aa = ^coth(V7WJ ( 35 ) 
i=i 

M a R = — — - if (a/3) is a bond (36) 
= otherwise , 
we may express the determinant as 

s (7) - _ n 2^u> I ( n dM*.)*-*=-»^ . (37) 

After performing the integration over </> Q one finds 

5(7) = II sh (V7W) det ( M ) (38) 

which is the expression given in |T]J (see appendix [A]) up to an inessential numerical factor 
(2ir) B ~ v that will be dropped in the following. 



5 The spectral determinant in terms of the link matrix 

Equation (^) expresses the spectral determinant as the determinant of a finite V x V matrix. 
The matrix M describes the topology of the graph, telling us which vertices are coupled, and 
also contains the metric information. At this stage, the metric and topological information are 
inextricably mixed. The purpose of this section is to derive another expression of the spectral 
determinant in terms of the determinant of a 2B x 2B matrix whose natural basis is the set 
of arcs (oriented bonds). An advantage of this formulation is that it leads to an expansion of 
the spectral determinant as a sum of terms that will be interpreted in the next section as the 
contribution of periodic orbits. 

Our starting point is expression fl34|). Since the 7 dependence of the determinant is sim- 
ple, one may set 7 = 1 and recover the 7 dependence at the end. The first step is to find 
a more convenient expression of the propagator in (B3). To begin let us consider the one- 
dimensional harmonic oscillator propagator gt(x,x') = (x\e~^~ dx+x ^\x'). This propagator 
may be expanded over the eigenstates <p n (x) = 1/4 ^ r -^ ff n (x)e~2 :i:2 of energies E n = n + |, 
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where H n (x) are Hermite polynomials. Using the generating function of Hermite polynomials 
S^=o H-n(x) = e 2Xx ~ x , the propagator may be rewritten as 

00 n -nl 



'a/3 



91 



a/3 



n=0 

A re-summation of the series gives 



v ' p; ^ 2™ re ^ 



A a /3 



(39) 



7T 



(40) 



Writing the two-dimensional propagator as a product of one-dimensional propagators Gi a „ (f a , rp) 
gi a0 {x a ,x/3)gi al3 (y a ,yi3) one is led to an analogous expression for Gi ag (f a ,rp). For each bond 
one has to introduce two couples of variables (\ a p, Xp a ) and (A^, X'g a ) associated to x a , xp and 
VaiVp respectively. Using the complex notations z a = x a + \y a and k a p = A Q/ g + iA^g (one 
recalls that d z = \{d x — id y )) one eventually finds 



7T 



where O a n is the following operator: 



O a(i = exp (e- l ^(d AafJ d Apa + d Aa3 d Af J 



1/3 



(41) 



(42) 



Equation ( |4l"| ) shows that one has to introduce two variables K a p and hp a per propagator in 
(|34]), i.e. each arc (a/3) is associated with a A Q/ g. 



From (|34|, 41 ) it follows that 



S( 1 = l)~ 1 = 2 



B-V' 



O 



7T* 



(a/3) 



(43) 



where the operator O = II (0/3) ®a/3 is understood as acting on both terms that follow it (one 
has multiplied ([34|) by the numerical factor (2ir) v ~ B to make it disappear in the final result). 
Again, one has replaced in (|J|) the product over bonds by a product over vertices. Integration 
over (f) a is then straightforward and leads to 



5( 7 = 1) 



-1 



2^e 



Oe 



A+QA 



(44) 



where A is the 21?-component complex vector: 



A 




(45) 
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and 



At = ( A lai 



A 



A. 



2f3„ 



(46) 



in which the A a/ g are grouped by vertices. The first group (Ai Ql , • • • , Ai a ) is related to the 
mi arcs ai, c*2, • • • , a mi issuing from vertex a = 1. In this basis, Q is a 2B x 2B block diagonal 
matrix (with V blocks); it only couples arcs that start from the same vertex 



Q = 



( Qi 










o \ 





Q2 
















Q 3 

















V o 










Qv J 



(47) 



where each m a x m a sub-matrix Q a attached to a given vertex a has ^ 1 on its diagonal 

and — everywhere else. Using a tensorial notation, one may write the matrix element of Q 
between the two arcs (a(3) and (fiv) as 



Q(a/3)(pu) — a a pa^ u 5 afl ( 5p v 



(48) 



where 8 a p is the Kronecker symbol; the connectivity matrices a a/ 3 and a^ v are here to recall that 
(a/3) and ([iv) are arcs. 

The action of O on e At « A is complicated because the argument of the exponential is quadratic 
in the A Q/ g's. This action would become much more simple on the exponential of a linear form; 
this can be achieved by introducing an additional integration over a 2£>-component complex 
vector W 



5( 7 = 1) 



2 B e~ L 



(Ha 171 



- T ^det(Q- 1 ) / dW^dWe~ wt ^ lw Oe Atw+w ' A . (49) 

a) ft J 

The action of O is now simple by noticing that d Kafj d K ^e k]w+w] K = w a pWp a e^ w+w ^ K . To 

write the action of O on e Atw/+w/tA , one has to replace d^ afj by w a p and by w a p. One 
finally ends up with 



r 2 * e * detCQ" 1 ) [dW^dWe- w ^- lw+w ' RW , 



^(7 = ir 1 = 



(EL™ 



(50) 



where it is clear from the expression of O that the matrix R only couples the arc (a/3) to the 
time reversed arc ((3a). Thus, R may be written 



R 



(51) 



It is then straightforward to perform the Gaussian integration and get the final expression for 
the spectral determinant 



5( 7 ) = ^e^ L {j{m a )2' B det(l - E) . 



(52) 



where we have introduced E = (QR) T (the transposition is defined as (-E? T )( a/ 3)(^) = ^(^(a/?))- 
Using this tensorial notation the product QR is particularly simple to perform 



(53) 
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The spectral determinant is now given by the determinant of a IB x 2B matrix E whose basis 
consists of the set of arcs. Rewriting the matrix E = De one may clearly separate a part e that 
only depends on the topology of the graph 

e (a/3)(pu) = a a pa^ u 5i3^ ( 5 av J (54) 

and a diagonal part D depending on the metric properties 

D {af3)(pu) = a a f3a^a^l3uQ~^ l>1V ■ (55) 

We end this section by mentioning the physical meaning of the different matrices that appear 
in the derivation given above. 

In order to simplify slightly our notations arcs will now be labeled by Roman letters . . . 
(i = (a/3)). The matrix E is then expressed as Eij = eije~^ li . Our ey coincides with the one 
introduced in [21]: = ^- if the end of arc i and the beginning of arc j are the vertex a of 
coordinence m a , ey = ^ 1 if moreover the two arcs are time-reversed, and = otherwise. 



Following Kottos and Smilansky [TJ|], one may take a scattering point of view to construct 
the wave function on the graph. Let us consider first the set of arcs i G {1, ...,m a } starting at 
vertex a. On each arc i one writes the wave function as the superposition of an incoming and 
an outcoming plane wave (see figure ||): 

^(i) = %- ifa + Bie ifa . (56) 

The scattering matrix s a at the vertex is defined as the matrix relating incoming to outcoming 




Figure 2: 

amplitudes: Bi = J2j( s a)-ijAj for i,j G {1, ...,m a }. Conditions ( p6| , p7| ) lead to the following 
result: s a = Q a . On the graph, the 2B outcoming amplitudes are then coupled to the IB 
incoming amplitudes by [13] 

B i = ^2Q ij A j (57) 

3 

(matrix Q indeed couples arcs beginning at same vertex). The physical meaning of the matrix 
Q is now clear: it describes the scattering by the V vertices. As in section || we use the notation 
i for the time-reversed arc of i. Amplitudes of two time-reversed arcs are related by the two 
obvious relations: Bi = Aie~ ikli and Ai = B^e lkl1 (due to the fact that ipi(x) = —x)). Then 
one has Bi = ■ Qije lkl ^Bj. Using the fact that Q-^ = one finds that 

B- i = J2E ij (^ = -k 2 )B- j . (58) 
j 

Then the state (|56|) of energy k 2 belongs to the spectrum if det(l — E{^ = —k 2 )) = [20|. As 



equation det(M(7 = —k 2 )) = 0, this equation gives the spectrum of the Laplacian but not the 
7-dependent pre- factor of 5(7) . 
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6 A trace formula 



The result ( p2l) is particularly suitable to derive a trace formula that expresses the determinant 
in terms of contributions of closed trajectories (orbits) on the graph. Again we may set 7 = 1 
for the sake of simplicity. We may expand the determinant 

00 1 

lndet(l-E) = -^-Tr{£ n } (59) 

n=l n 

and write the trace as a sum of terms 

Tv{E n } = e ilia e iai3 ---e inil e-toi+- +l *\ (60) 

each of which is associated with a sequence of arcs that forms a closed path on the graph. 
Following the notations of [21] one may attach the quantity 

Ot{Cn) = ■ ■ ■ i.i n i x (61) 

to the orbit (or circuit) C n = (ii,i2, • • • >&n) composed of n arcs, and denote its length by 
l(C n ) = li 1 + ■ ■ ■ + li n . Each term of TV {E n } corresponds to a path associated with either 
a primitive orbit of n steps or k repetitions of a primitive orbit of p steps such as n = hp. 
Primitive orbits will be labeled as C. A given primitive orbit of n steps appears n times in 
Ti{E n }, corresponding to the n cyclic permutations of the indices that give the same orbit 
on the graph. On the other hand the term associated to k repetitions of a primitive orbit C p 
appears p times corresponding to the p different permutations of the indices that appear in 
Tr{E n }. Using these remarks, instead of summing over all paths of length n as in (|60|), one 
may sum over the fc-repetition of primitive orbits of length p provided that n = kp: 

Ti{E n }= J2 P a(C p ) k e- kl ^ , (62) 

k,C p with n=kp 

we have used the fact that if C n is a ^-repetition of C p then ct(C n ) = a(C p ) k and l(C n ) = k l(C p ). 
Introducing this expression in (59) gives 



00 ^ 

lndet(l-£) - ^Er^^^' ' ( 63 ) 
c P *=i P 

Instead of summing over n, k and C p with the constraint kp = n, one may sum over all primitive 
orbits C p and their /c-repetition independently. This immediately leads to 

In det(l -E) = J2 ]n ( 1 ~ a(C p )e- l( -^A . (64) 



One thus gets an expression of the spectral determinant in terms of an infinite product over all 
the primitive orbits C, in infinite numbern, that can be constructed on the graph: 



5( 7 ) = ^e^ L mm a )2- B l[(l - a{C)e~^ 1 ^) . (65) 



c 



6 Note however there are two trivial graphs for which one can construct only one primitive orbit: the circle and 



the line. 
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A similar trace formula for the partition function was first derived by Roth [21]. It is indeed 



possible to recover it by starting from (65). The log-derivative of the spectral determinant is 



equal to the Laplace transform of the partition function (0). From (65) one gets 



fr™-£ + I ir + £jX>®± ««5)'e-^<*> , (06) 

where one has expanded the denominator coming from the derivative of the logarithm that 
appears writing In £(7). It is then possible to group the two sums in a unique sum over all 
orbits C, including the repetitions of the primitive orbits 

#- In S( 7 ) = J- + Y —^- + J- £ l{C)a{C)e-^0 . (67) 

In this summation, C designates the primitive orbit associated with a given orbit C. 

p. 

Using the following identity ^=e~^ 1 = -^j^ / °° dt -^e _7 *~4t it is then easy to extract the 
inverse Laplace transform of the previous expression and recover the trace formula first obtained 



by Roth [21, 22 



zlt) -ik +v ^ + ^sV e)aiC)e ~^- (68) 

Let us notice that if one considers a vertex a of coordinence m a = 2, the circuits that 
contain a reflection at this vertex have weights a(C) = 0, therefore such circuits should not be 
considered in the expansion. Since one can always introduce an arbitrary number of vertices on 
any bond without changing the properties of the graph, this remark insures that the number 
of orbits will not vary doing so. Moreover this shows that a graph can always be simplified to 
minimize B and V by suppressing all vertices of coordinence 2. 



7 Simplification of the trace formula - Diagrammatic Expansion 
of the spectral determinant 

In the previous section one has expressed the determinant det(l — E), which can be expanded to 
give a finite number of terms, as an infinite product (|65|). The purpose of this section is to show 
how the infinite product eventually simplifies to give a finite number of terms. One interest of 
this discussion is to provide a diagrammatic method to construct systematically the different 
terms of S(-y). 

We will consider the quantity that appears in ([65]) : 

S = Yl(l-a(C)e- l{d) ^ . (69) 

c 

The generalization of the following discussion to the case with a magnetic field is straightforward 
(see ©). 

We will call — a(C)e~ 1 ^ the weight of the orbit. We associate to each term in the expansion 
of S a diagram, that represents either the contribution of an orbit or the product of such 
contributions; in the last case the diagram represents the superposition of the different orbits. 

To understand how the simplifications occur let us consider a primitive orbit C that passes 
twice through the same vertex without using the same arcs as on figure [|[ Using the fact that 
£12 = ^34 = ^14 = £32 it is easy to see that 

_ a(C)e- l & = - (-a(C")e-^')) (-a(C")e-^")) (70) 
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Figure 3: 

where & and C" are obtained by crossing the paths at the vertex as shown in figure ||. The 
consequence is that in the expansion of S, the product of the weights of orbits C and C" cancels 
with the contribution of C. 

These kind of relations may be represented diagrammatically as 




Using these relations one must take care of not introducing some reflection at a vertex like in 
figure |j. Indeed the diagram on the left of figure || has a weight proportional to (^-) 2 whereas 
the diagram on the right has a weight proportional to — (— 1) — • 




Figure 4: Take care not to introduce some reflection at a vertex using the rules (|7T|). 

Note that the case of a vertex with coordinence m a = 4 is special and may bring some 
additional rules which simplify the expansion of S. Those two additional rules are: 



and 




What kind of consequences can be deduced from these rules ( |7l| ) ? Let us consider a primitive 
orbit that contains twice the same arc. This orbit may be factorised using the de-crossing rule 
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(74) 



This equation implies that one should not include in the expansion of S the diagrams that 
contain more than once a given arc. This remark implies that the number of diagrams to be 
considered is finite. In the diagrams that remain in the expansion, a bond appears in the orbits 
at most twice, corresponding to the two reversed arcs; a by-product of this remark is that the 
longest orbits that can be constructed satisfying this rule have lengths 2L. 

Example of a diagrammatic expansion 




2 



Figure 5: 

As an example, consider the diagram of figure [|, which consists of an arm connected to a 
ring pierced by a flux (ft; this geometry was considered for the study of persistent currents in the 
case of a one-channel clean ring |33| and in the case of a metallic diffusive ring p| . The first 
step is to construct all the periodic orbits that will be involved in the expansion (see the table 
|lj). One considers the situation where the ring is pierced by a flux <f> to distinguish between 
time-reversed orbits; = 2~K<f)/<j)Q where cfto = h/e is the flux quantum. We call I the perimeter 
of the ring and b the length of the arm. As an example let us compute the a(C) coefficient for 
the last orbit of the table |l| C = (1, 2, 2, 1): a(C) = ei2e22 e 2i e ii with en = — |, eu = egj = § 
and e 2 2 = 1- 

The construction of all diagrams of the expansion requires to combine all these orbits provided 
that the resulting diagrams do not contain twice the same arc. This leads to 




where it is understood that a diagram with two primitive orbits is equivalent to the product of 
diagrams for each orbit, as said above. For example 



(76) 




Using the weights of the orbits given in the table one eventually finds 



I — - cos ( 



cos ( 



-l-2b _j_ e — 2Z— 26 



(77) 
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primitive orbit C 



weight —a(C)e 



-l{C)+i6{C) 




2-1-18 
"3 e 



_g -i+ie 

3^ 



1 Q -2Z 



4 a -l-2b-i8 



4 -l-2b+W 
'9 e 



4_„-2l-2b 
27 e 



J^„-2Z-2fc 
27 e 



Table 1: 



-l-b 



(sh 6 sh Z + 2(ch I - cos 6) ch b) 



Thus: 



5( 7 ) = sh(V^>) sh(V7i) + 2[ch(V70 " cos(0)] ch(^>) , 
a result that can also be obtained using ( |35| , |38| ) (see also appendix 0). 



(78) 



(79) 



As one can realize looking at ([77D, there exists a symmetry between the coefficients of the 
different terms appearing in the diagrammatic expansion. It is possible to use this symmetry to 
reduce the number of diagrams that have to be considered, which greatly simplify the calculation. 
This point is discussed in more details in appendix [B|. 

To conclude this section we insist on the fact that the periodic orbit expansion of the spectral 
determinant involves only a finite number of contributions, despite the number of primitive orbits 
is infinite. In contrast, the periodic orbit expansion of the partition function (|l0|) or the density 
of states involves an infinite number of contributions. 



8 Graphs in a magnetic field 

In this section, we describe the appropriate modifications to the previous formalism that have to 
be done in the presence of a magnetic field, namely for a distribution of Aharonov-Bohm fluxes. 
The operator of interest is now — (d x — L4) 2 . In equation ( p7[ ) one has to replace the derivative 
of the function tp by a covariant derivative D x = d x — iA(x) 

^2 a a pVx al3 i>(aP) (Xa/3 = 0) = . (80) 
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The path integral derivation of det(— + 7) follows the same lines except that all derivatives 
have to be replaced by covariant derivatives in ([30],[n]). One is led to an expression of the 
determinant that involves a product of terms of the form 



0(O)=0< 



We4^ dl fl D ^ 2+ ^. (81) 



We may recover the Gaussian action of the harmonic oscillator by performing the following 
gauge transformation 

<t>{x) =4>(x)e iJ *o dx ' Aix ' ) (82) 

where the integral is performed along the bond (a/3) and xq is an arbitrary point on this bond 
(a change of xq corresponds to adding a constant phase to the field). If we define 



'Pa 



^ dxA(x) (83) 



(Qpa = ~®ap) where the integral is performed along the bond (a/3), and choose xq in such a way 
that fj*> dxA(x) = % and f° Q dx A(x) = then @ becomes 

r4>(ic,0)=<f>0e~ i<>lSa/2 -, 1 a 

fl #e4/o a ^(MW) = ^ G (7 i/4 0/3e -i^/2 )7 i/4^ e -i^/ 2) _ 

(84) 

Introducing this expression in d34|), then it is clear that (^) still holds provided the matrix M 
is now defined as 

M a/3 = 8 a p °V coth( v /7/ Q , At ) - a a/3 — r . (85) 

sJl l\/7'a/3j 

The calculations of section || may also be generalized. The derivation is the same and the 
matrix Q now reads 



Q(ptfi)(jw) — a a(3 a ^an ( 5p v j e ' 2 ! (86) 



and the matrix i? is unchanged. To define E 1 it is convenient to introduce a unitary transforma- 
tion U that makes the changes more clear. If one writes 

QR = UE T tf, (87) 

where U is 



then one has -E = De where e is still given by (|54| ) and all the dependence in the fluxes is now 
contained in the matrix D: 

D (a/3)(nv) = a aP a^5 a ^p v e-^ l ^ + ' e ^ . (89) 

Since the spectral determinant is given by the determinant of the matrix 1 — E, the matrix U 
and its inverse disappear in the determinant and ( |52"| ) still holds. 

The derivation of the trace formulae only requires trivial modifications and one finally obtains 
for the determinant: 

5( 7 )= 7 ^ e ^ i (n^)2- B n( i - a ^) e ~^ /(c5)+i ^ ) ) > ( 9 °) 

a q 
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and the corresponding partition function: 

m = i + ^ + ncMOe-^v (91) 



2V^i 2 2V^i c 



which generalizes the Roth's formula (|68|). We have used the obvious notation for the flux 
enclosed by an orbit C = i%, ■ ■ ■ , i n ) : e(c) = e h + ■ ■ ■ + e in . 



9 More general boundary conditions 

One has sometimes to consider more general boundary conditions than (|30|). Imposing instead 
of ( |80|) the so-called mixed boundary conditions 

^2 a a f^>x a ^{ a p) {Xaf3 = 0) = \ a Tp a (92) 

p 

requires to generalize some of the previous results]]. This can be easily achieved by using the 
path integral formalism. It is rather obvious that the boundary terms in (^Tj) now produce 
additional Gaussian term in which now reads 

r v a 

S( 7 )-i = 7^ / II d^d^e"^ 1 ^ 12 J] G^^e^W^e-^) . (93) 

J a=l (a/3) 

Thus only the diagonal elements of matrix M change. In (]38|), M is replaced by M: 



Map = M a/3 + ^8 aP . (94) 

It is also easy to see from the discussion of section |5| how the different matrices are affected by 
the change of boundary conditions. Equation ( |43| ) will receive some additional Gaussian terms 
e - 2 Aal^al 2 /^ whjch implies that m a has to be replaced by m a + The scattering matrix Q 
now reads 



and the matrix e is 



e (a/3)(pu) = aa/3<V<W( TT ~ ^auj , (96) 

the other matrices R and D remain unchanged. The factor a(C) of sections |6] and |7| will be 
modified according to (|96|). 



Expression (94) allows us to treat both the case of Neumann boundary conditions (A Q = 0) 
which has already been studied in this paper, and the case of Dirichlet boundary condition 
(Aq, = oo). For a diffusive conductor, these conditions correspond to a disconnected wire and a 
wire perfectly connected to leads at the node, respectively. 



7 One may distinguish two cases, (i) In the Schrodinger problem, the current is J = 2Im(ip*d x tp). Conditions 
( pif ) also lead to current conservation at the vertices, (ii) For the diffusion problem, the wave function ip(x, t) has to 
be replaced by a probability density P(x, t). The current of probability is given in this case by J(x, t) = —d x P(x,t). 
Equation (B2|) implies that such a current is not conserved at the vertices where A Q 7^ 0. This condition corresponds 
physically to a graph connected to the external world, \ a P(a, t) being the current of particles exiting the graph. 
In the Dirichlet case (X a = 00) the graph is perfectly connected, which means that a particle exits the graph with 
probability 1 if it reaches the vertex a. 
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If one imposes Dirichlet condition at all vertices of the graph, one has to consider the 
limit A — > oo for all vertices. Then det(M) ~ Y\ a The spectral determinant is S(j) ~ 

( Yl a 11(0/3) sh ^ Q<3 • One can drop the irrelevant factor Y\ a X a since the spectral determi- 
nant is defined up to a multiplicative numerical factor independent of 7, that depends on the 
regularization. One finds that the spectral determinant 5(7) = E[(a/3) sh ^ a/3 is the product of 
the spectral determinants associated with each bonds. All bonds are then independent. 

We now discuss the case where one imposes Dirichlet boundary at only one vertex ao which 
will be useful in the next section. Taking the limit X ao — > 00, one has: 

5(7; {Dirichlet at q }) = 7^^ JJ sh(^ a/3 ) det(M Q °) (97) 

{a/3) 

where M a ° is the (V — 1) X (V — 1) matrix given by the matrix M with the line ao and the 
column ao deleted. 

As an example consider the graph pictured in figure ||| for which one imposes Dirichlet 
boundary at the end of the arm. S(-y) becomes: 

S Dirich . at (7) = 4={ ch (v^ & ) + 2[ch(^70 - cos(0)] sh(^7&)| , (98) 

end of arm y 7 ^ 



to be compared with (79). It is interesting to check that the magnetization of the ring, given by 



(1^,2^), does not depend on the choice of the boundary conditions at the end of the arm, when 
its length b goes to infinity (b 3> = 1/ y^y) . 

10 Zero mode and low energy behaviour 

The low energy part of the spectrum can be studied by using the expansion of the spectral 
determinant when 7^0. Using this approach we prove the existence of a zero mode state 
when B = and obtain the dependence of the ground state energy in the magnetic field. Using 
the infinite product representation SM = Yl n (En + 7) an< ^ assuming Eq = implies that S('y) 
behaves linearly with 7 when 7^0. In order to study this limit it is convenient to expand M 
as a power series in 7 

MM = — M(7) = — (M° + 7M 1 + ---) (99) 



where 



= i^Er-r (100) 



'a/3 

Ms/3 = ^3 a "/^"M + g a a/3^a/3 • (101) 
I" 

From these expressions, it is easy to see that M° possesses an eigenvalue £o = corresponding to 
the eigenvector vq whose components are (fo)a — 1- Remembering that M acts on the V^-vector 
constructed with the wave function at the nodes, the vector vq corresponds to a wave function 
that takes the same value on all vertices and is constant on the graph (since 7 = —E = 0). We 
denote by ^(7) (n = 0, • • • , V — 1) the V eigenvalues of MM. The eigenvalue £0(7) vanishes at 

7 = and its behaviour for small 7 computed in perturbation theory is £0(7) — 7 "° ^ V ° = -rr- 

Vq Vq v 
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This implies that det(M(7)) ~ ^ Iln=i £n(0)- It is possible to write the product of the non- 

T 

vanishing eigenvalues for 7 = as the determinant of the matrix M° + where the second 
term is the projector on the vector vq. 

One finally obtains for the spectral determinant 

S{i) ^^(Hl aP )det(K) (102) 



7- 



where the matrix K is 



1 



K a/3 = M^ + -. (103) 

From the linear behaviour £(7) oc 7 at small 7 it follows that one has for the partition 
function: lim^oo Z(t) = 1. 

Let us remark that if all the lengths of the bonds are equal to unity, the product n„=i £n(0) 
is equal to T(Q), the number of trees that cover the graph (Tutte theorem) |l6], [T^] . 

It is also interesting to see how the spectral determinant behaves at non vanishing magnetic 
field. For the sake of simplicity we denote the ensemble of fluxes {0 a p} by B. The matrix M now 
depends on the fluxes. The spectral determinant may be expressed in terms of the eigenvalues 
of M(7,B): 

S(l,B)=l\ S -^^liui,B). (104) 

(a/3) n=0 

At zero magnetic field and small 7, one has shown above how S("f) behaves ( |102| ) by a pertur- 
bative expansion of Co (7)0)- Similarly, at small magnetic field, one may compute the eigenvalue 
£0 (0, £>) in perturbation theory. Starting from 

M a/3 (0, B) = 5 a(3 J2r t - a -/3^ ( 105 ) 

,, 'a/z 'a/3 

02 

and expanding this expression at small fluxes, it is easy to show that £o(0, B) ~ ^7 J2 a p a Q /3-p^- 
Since 11^=1 Cn(0,0) = det(K) the spectral determinant behaves like 

s(o,{o«p}) * iEr(IIw) det W' ( 106 ) 

One may extract from these expressions the ground state energy at small magnetic field. At 
small 7 one has: 5(7,0) ~ 7 Yl n>0 E n (0) . On the other hand, when the magnetic field goes to 
zero, the ground state energy is the only energy to vanish, then S(0,B) ~ Eo(B)Y\ n>0 E n (0). 

The ratio — gives the expression of the ground state energy 

Eo (R/3}) ^Er' ( 107 ) 



Let us notice that the expression (107) of the ground state energy at small magnetic field 



may be recovered in a simpler way by perturbation theory. The eigenfunction of the ground state 
at zero field is tpo(x G Q) = -7=. Thus the correction to the energy at small field is Eo({9 a p}) ~ 

( V>o I^K^) 2 ! i>o ) = 7j Yl( a /3) lo™ 13 ^ x A(x) 2 . By choosing a gauge such that A(x) = Ap a is constant 
on the bond (a/3), one recovers ( |107[ ) since 9p a = Ap a l a p. 
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So far we did not describe the effect of a change of boundary conditions on the low energy 
behaviour of the spectral determinant. The matrix M at zero magnetic field now becomes 



M a(3 (0, {X a }) = M% + X a 6 aP . (108) 

As in the presence of a small magnetic field one may compute the lowest eigenvalue £o(0, {A a }) 
of the matrix M(0, {A Q }) by perturbation theory. This gives 

S(0, {X a }) ^ i V A J 17 Z a/3 ) det(iO • (109) 
One may also deduce from this result the behaviour of the ground state energy: 



L 

a 



Let us notice that in a particular case, one may recover this result by another way. Consider 
a ring with V vertices with mixed boundary conditions. This problem is equivalent to a ring 
with a potential made of 5 scatterers: W{x) = Yl a ^a5(x — x a ), where the the positions 

of the vertices on the ring. In this case it is possible to compute perturbatively the energy of 



the ground state: Eq {{Xa}) ~ (i/)q \W(x)\^q) which indeed leads to (|110|) . 

The result ( |110|) has a physical interpretation in the context of diffusion. The zero mode 
Eq = is associated with the uniform stationary distribution at long time. As soon as the graph 
is connected to the external world, which may be described by taking mixed boundary conditions, 
the probability to remain on the graph decreases at long time J" Gra h dx P(x, x' , t) ~ Q~ E ot which 
means absence of a zero mode. In this case, Eq is the inverse escape time of the graph. 



11 The graph connected to an infinite lead 

In this section, we consider the case of a graph connected to an infinite lead attached to the 
vertex a. In that case, we are dealing with a quantum scattering problem where the relevant 
information like the spectrum, the Wigner time delays or thermodynamic quantities like the 
persistent currents are encoded in the (unitary) scattering matrix |34|] . To the purpose of 
calculating the scattering phase shifts, let us suppose that along the lead an incoming plane 
wave e~ lkx enters the graph at the vertex a and is reflected with a phase shift 5 a (E). Below 
we derive a formula expressing the phase shift in terms of det(M). Although such a formula 
is in fact a special case of the one given in [l4|], it is nevertheless interesting to present 
an independent derivation of it and show the relation with spectral determinants. m a is the 




Figure 6: 



coordinence of vertex a in the absence of the external lead. We first write the wave function of 
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energy E = k 2 on the bond (n/3) in a form that insures continuity: 

(x^fi) = . U1 4>v sin k{l a p -x^)+ tppe lA ^ 1 ^ sin kx^ (111) 

sin K>(>ij J i3 ^ * 

where AruIu/3 = On the incoming lead, the wave function is written in terms of stationary 
scattering states 

Vw(*) = -4fl^ cos ( kx + <W 2 ) « e~ ikx + e ikx+iS « (112) 
cos(o a /2J 

where x G [0, oo[ is a current point on the lead. Continuity of the wave function at vertex a is 
already insured with the previous expression. Current conservation at the vertex a gives 

2 , a a/3 I COtg /c/ a/3 + i>p-^-j-j ~ Ipa tSn(6 a /2) = . (113) 

Since the boundary conditions on the other vertices are obviously unchanged we are left with 
the linear system of V equations: 

- i Mtfb = -k 2 )ip p = 5 ua ip a tan(<5 a /2) (114) 

where fi = 1, V. Cramer's formula then gives 

fe=itmfe /2) jg^» fe (115) 
where M a (7) is the matrix introduced in section ^. Then the phase shift is obviously given by 

This result, together with (|97|), shows that the phase shift is related to the ratio of two spectral 
determinants: one calculated with a Dirichlet boundary at the vertex connected to the lead and 
Neumann boundary at all other vertices (|97|), and the other calculated with Neumann boundary 
at all vertices (Q): 

cot g( ^)/2) = . (117) 

This expression shows that 5 a = mod 2ir if k 2 coincides with an energy of the graph with 
Neumann boundary, and 5 a = it mod 2-k if k 2 is an energy of the graph with Dirichlet at vertex 
a. 

As an example the phase shift of the graph of figure ||| is given by the ratio of © and © 
if the lead is connected to the end of the arm 

s cos A:6sin kl + 2(cos kl — cos 8) sin kb 
CO,B(i/2) " 4n<**n<U -MtU,* ' (U8) 

As a second example the phase shift in the case of the complete graph K n that will be studied 
in the next section reads 

ootg(*/2) = cos( V? ) COS( ^ + 7^. ) " 1 cotg^/2) , (119) 
cos(ki) + cos{ip) 

where co$((p) = —^r. As expected the scattering length 4f | fc=0 = 2 n ^ w ~ 1 - > l = 2L is a measure of 
the total length of the graph. 
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Figure 7: The complete graph K§. 



12 Application to the case of the complete graph K n 

As an illustration of the previous formalism it is interesting to consider the particular case where 
the graph Q is the complete graph K n whose V = n vertices are all connected (see figure [?]). 
Then the number of bonds is B = n~ ■ 

In the case where all the bonds of a graph have the same length I one obtains 



det(M) 



1 



P((n -l)ch(V^)) 



(120) 



where P(X) = det(X5 a /3 — a a p) is the characteristic polynomial of the graph. The characteristic 
polynomial of K n can be found in classical textbooks. For the complete graphs a a p = 1 — 8 a g, 
thus the determinant involved in P(X) is of the same form as the determinant of matrix F given 
in appendix |I^. This gives the following expression for the spectral determinant 



S{l) = {Vl)~~~ (sh Jll)-~-(n- l)(ch ^y£ - 1) [(n - 1) ch ^ji + 1] 

d_ 

d-, 



n-l 



(121) 



By calculating the inverse Laplace transform of -S-hxSipf) one obtains the partition function 



_ n(n-l)£ w(3 - n) 
AV^i 4 
^n(n — 3) 



2V7rf 
wt 



V 4vrt 



1 



+ 2 (n - 1) £(-l) fe T fc 



fc=i 



n — 1 



e 4t 



122) 



where 9(x) = Y^t=-oo ex P ~ 7r k 2 x is the Jacobi function, and T^{x) Tchebychev polynomials. 



a) short time limit 



Equation (|122j ) is written in the form of Roth formula, appropriate for short times. The first 
contributions coming from the expansion of (|122j) in the limit t — > give 



Z(t) 



L V-B 1 
+ — = — + 



2Vvrt 



n(n — 1) 
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1 



n 



e t 



n(n — l)(n — 2) „ 
+2^ ^ ^3£ 



n 



e « + 



(123) 



The first two terms coincide with the first terms of Roth formula. For the two following terms 
we recognize the contributions of the shortest orbits represented in figure M. For the orbits of 
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length 2£, the factor — lj is the weight a{C) of those orbits and " ^ l ' their number. 

The next term corresponds to orbits of length 3£ with weight a(C) = (^zj) • Their number is 
2 «(n-iKn-2) where the additional factor 2 comes from time reversed orbits. 




Figure 8: Some orbits associated with the first terms in the brackets of ( |123| ) for the complete 
graph K 6 . 



b) Long time limit 



For the graph K n it is also possible to study the limit t — * oo. Introducing cp, defined by 
cos ip = — , one may use the identity Tk (cos (p) = cos kip and the Poisson summation formula 
to find 



+ {n-l) e 



(124) 



fc=— oo 



It is easy to see that 

Z(t) = 1 + O (e~ Kt ) (125) 
where the contribution of the first term comes from the constant zero mode ipo(x £ Q) = -^=. 

c) The eigenvalue spectrum 



Equation (124) provides the whole spectrum of K n . The first two terms describe a series of 



states E\2,k = k' z jr, k £ N with degeneracies c£i2,fc=o = 1> d\2,k 
d\2,k = 2 + "^" 2 3 ^ if k is even. The third term generates energies Ez,k = jr (2& — 1 — f ) 2 , k £ Z 
with degeneracies tfg^ = n — 1. These states are similar to those of a ring pierced by a flux, 

where the parameter (p related to the coordinence plays the role of the flux. 

rill 



n(n— 3) 



if k is odd and 



For n ^ 3, ^ £ [|, |[ and one may reorganize the energy levels: 



Enerj 


;y 






Degeneracy 


So = 


Sl2,0 = 


= 




rfo = l 


Si = 


^3,1 = 


7T 2 

w 




d\ = n — 1 


E2 = 


Si2,l = 


7T 2 
= W 




j n(n— 3) 


S3 = 


^3,0 = 


7T 2 

~¥ 




c?3 = n — 1 


S 4 = 


Sl2,2 = 


7T 2 

= "F 


2 2 


d 4 = 2 + n(n - 3) 


s 5 = 


^3,2 = 


7T 2 
T 2 " 


(3-f) 2 


ds = n — 1 
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The case n = 3 corresponds to the ring of length 31. In this case ip/ir = 1/3 and one recovers 
the well-known spectrum. 

It is interesting to note that despite the particle explores a volume L ~ \n 2 l at large n, the 
energy of the first excited state E\ ~ |w is not of order j?, as one could have naively guessed, 
but instead of order -p. 

13 Conclusion 

We have investigated spectral properties of the Laplacian on graphs by providing several equiv- 
alent representations of the corresponding spectral determinant. This has been achieved thanks 
to a path integral formulation and to the fact that the spectral determinant can be written in 
terms of the propagator of a 2D harmonic oscillator thus leading to a set of straightforward 
Gaussian integrals. We have thus obtained an expression of the spectral determinant as a trace 
formula involving the contribution of an infinite number of periodic orbits. Using a systematic 
diagrammatic method, we expressed the finite number of terms of the determinant as contribu- 
tions of a finite number of orbits. Although it has been already clear in the literature that such 
a reduction is possible, the present formalism allows us to implement it directly for any given 
graph. 

The flexibility and the relative simplicity of this formalism is a hint for using it in a broader 
range of problems. For networks of mesoscopic and coherent conductors, it may help to compute 
both thermodynamic and transport properties in some non trivial situations such as the local 
distribution of persistent currents in a network driven by a far remote Aharonov-Bohm flux. 
Along the same lines, it could be possible to address the problem of topological properties of 
fractal networks like the Sierpinski gasket subject to a distribution of Aharonov-Bohm fluxes. 
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A The spectral determinant obtained by constructing the Green's 
function on the graph 

We want to solve the diffusion equation 

{ 1 -&l)G(x,y)=5(x-y) (126) 

on a graph made of V vertices (or nodes) linked by B bonds, y is the source for the diffusion. 
The solution of ( 126 ) is of the form: 

c(^) = £*^|M (127) 

where i\) n are the eigenfunctions of the operator — d^. The spectral determinant is obtained by 
spatial integration of the diagonal Green's function G(y, y) : 



dyG(y,y) = £ -J— = ^-lnS( 7 ) . (128) 



7 + E n 
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For a given source located in y, the diffusion equation is solved on each bond (a/3) in terms 
of the values G(a, y) and G(f5, y) at the nodes: 



G(x, y) = G(a, y) ch y^x + (G(P, y) - G(a, y) ch y/^l a p) 

where x is the linear coordinate on the bond (a/3) of length l a p. 
Current conservation at the vertex a: 



sh ^fx 
Sh <fyl a/ 3 



(129) 



~ Yl d x a/3 G(x a p = 0, y) = 8 a>y , 

13 

where the sum stands over the nearest vertices of a, leads to the following equations 

G(f3,y) _ 6 a , y 



G(a, y) Y coth ^ ~ Yj 
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(3 



shr/ Q/3 ^7 



(130) 



(131) 



where r] a p = yj^l a p. We have thus obtained a system of (V + 1) linear equations for the (V + 1) 
variables G(a, y) where a can be either a node of the graph or the source point y. 



( G(a u y) \ ( \ 



G(y,y) 



(132) 



\G(a v ,y) J \ J 



where the (V + 1) x (V + 1) matrix My is defined by equations (35|,36). Here, a and are either 
vertices of the graph or the source point located in y. One now wants to calculate G(y, y). First, 
using ( |131| ), it is written in terms of G(a,y) and G(b,y) where a and b are the vertices ending 
the bond to which y belongs: 



G(y, y) (coth rj ay + coth rj yb ) 



G(a,y) G(b,y) = 1 
shr] ay sh rjby ^7 



so that: 



(A ,1. 1 ) - ( Sh7?a C h ^ + G(a, y) sh Vby + G(6, y) sh ) . 



sh r] ab 



(133) 



(134) 



The variable G(y, y) can then be eliminated in the previous system. The two equations for 
G(a, y) and G(b, y) are modified as: 



G(a, y) coth ^ ~ Y G (P> y ^ 



11 sh rjyb 



sh rjap \fl sh T] ab 



(135) 



and the (V — 2) other equations are unchanged. We obtain now a system for the V variables 
G(a,y): 

I G( ai ,y) \ ( \ 



M 



G(a,y) 



G(b,y) 



\ G(ay_ 2 ,y) / 



y/ishr] ab 



sh r) yb 
shrjay 
V / 



(136) 
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where M is defined in (35,^6). a and b are now two of the V vertices and the source point y is 
now excluded. 

By inversion of the matrix M, one obtains G(a,y) and G(b,y): 



G(a,y) 



(T aa sh r] yb + T ab sh r] a y) 



(137) 



where T = M 1 . From (134), one finally obtains G(y,y): 

G(y, y) = ^—\ sh Voy sh + * — \rp aa gll 2 + sll 2 + 2 j- ba g^ ^ g^ ^1 I _ ( 138 ) 
yfi \ sh7] ab sh z r] ab L J 

The spatial integral of y) on the bond (afr) can be written as: 



dyG(y,y) 



^{??afeCOthT/ ab - 1 

+(T aa + T 66 )(l + 2 7 — ) cothr/ afe - 2T ba (l + ^-^-) — } (139) 



where we have used the equalities: 



•V 



sh 2 77 
chr; 
sn r/ 



27— cothn 
07 



2 7 - 



9 1 



9 7 shr; 

Summing over all the bonds of the graph and using the following identities: 



((Taa + T bb ) coth Vab - 2-^-) = Tr {TM} = V 

(ab) V ShT]abJ 



^ ( (T aa + T bb )^- coth rjaf, - 2T ba - 9 



(ab) 



3 7 shrj ab 



Tr <! T—M 

a 7 



5 

?7 coth rj = 2 7 — In sh n , 

07 

we find that the sum P = J Gra h dy G(y, y, ) simplifies considerably into: 

p = |^E lnsh ^ + ^ + Tr { M " 1 ^ M } • 

Using the following property Tr |m _1 ^;M| = lndet(M) one finally obt 



(140) 
(141) 

(142) 

(143) 
(144) 

(145) 



ams: 



P=^mS( 7 ) 



where 



5(7) = 7 " 



]Jshr7 afc det(M) 

(ab) 



(146) 
(147) 
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B The symmetry of the coefficients of the diagrammatic expan- 
sion 



The purpose of this appendix is to discuss the symmetry properties between the coefficients 
appearing in the expansion of the determinant discussed in section 0. This symmetry, a con- 
sequence of properties of the matrix E, considerably symplifies the calculation of the spectral 
determinant by the diagrammatic method. Since the presence of a magnetic field does not affect 
any point of the following discussion, one will forget it. One will start to consider first the case 
of Neumann boundary conditions and will indicate at the end of the appendix how to extend 
the result to the more general case. 

Let us write det(l — E) = det(-E) det(l — E^ 1 ). The inverse of E is easily calculated since 
e _1 = e T and D is diagonal: -D" 1 = 5ije^ li . Then one has 

det(l -E) = det(l - eD) = (-l) v - B e ~ 2 ^ L det(l - eD' 1 ) (148) 

(for the calculation of det(E'), see appendix The matrix el? -1 is equal to the matrix E in 
which one has replaced ypy by —ypy. This implies that det(l — eD) and det(l — eD~ l ) exhibit 
similar expansions (with the same numerical factors). 

To use this relation it is convenient to organize the expansion of det(l— E) in terms associated 
with diagrams or product of diagrams involving the same set of arcs. Let us call g n = {ii, ■ ■ ■ , i n } 
a set of n arcs of Q (in particular go = 4 an d giB = Q)- One may write 

2B 

det(l -E) = J2J2 K[9n]e-^ l[9n] , (149) 

n=0 g n 

where l[g n ] = In + • • • + k n - The coefficient n[g n ] is a sum of coefficients a{C) or a sum of 



product of such coefficients. Equation (148) implies that there is a relation between the terms 
involving diagrams with a large number of arcs and a small number of arcs. More precisely, as 
a consequence of the obvious relation l[Q — g n ] = 2L — l[g n ], we obtain 

k[G -g n ] = (-l) V - B K[g n ], (150) 

where Q — g n = {i/i ^ g n } is the set of all arcs of Q except those of g n . As a simple consequence, 
the term of det(l — E) related to the diagrams of length 2L is precisely (-l) v - B e~ 2 v^ L since 

The use of the relation ( |150| ) is particularly powerful since it allows one to consider only half 
of the expansion ( |149| ) . Moreover the terms K[g n ]e~^ l ^ gn ^ for small n not only involve a small 
number of diagrams but also the simplest diagrams, constructed with less than B arcs. Taking 
into account in only the B first terms in ( |149| ) considerably reduce the number of diagrams to 
be considered. 

As an example, one may apply this relation for the graph of figure [|. For the term n = 1 



in (149), the sum over gi brings four ensembles {1}, {1}, {2} and {2}. Only the two first give 
some contributions since {2} and {2} are not associated with any orbit. Consider the case where 
gi = {1}. Then the term K[{l}]e~v / 7 / [{ 1 }] i given by the diagram 

(151) 

has the same numerical factor as «[{1, 2, 2}]e~^ 1 ^ 1,2,2 ^ associated with the two following dia- 
grams 

(152) 
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In (|75| ) only the 3 first diagrams, among the 15, have to be considered and all other terms are 
given by the relation ( |150| ). It is also possible to check the symmetry ( |150| ) on equation (|77|). 

To end this appendix let us explain how the relation ( |150| ) is generalized for mixed boundary 
conditions. In this case the inverse of e is given by e T in which one has replaced ^7" by — y/7 
(the coefficients a(C) and n[g n ] now depend on ^7). One has det(l — E) = det(E) [det(l — 
e-D)|^y_^_^]. Then equation ( |150|) reads in the more general case 

k[G - g n ] = det(e) K\g n ]\^_^ , (153) 

det(e) being given by (|170| ). 



C Precisions on loops and multiple bonds between vertices 

In the section || we have assumed that a bond links always two different vertices and that there 
is only one bond between two different vertices. Formula ( |35| , p^ ) apply to those kind of graph. 
As we remarked in a footnote, if two vertices are linked by two bonds, one can introduce a vertex 
in one of those bonds to fall back to the situation of section || (see figure |9|). Similarly if a bond 
starts and finishes at the same vertex to form a loop it is possible to separate it into two bonds 
by introducing an additional vertex (see figure 10). Of course, these operations do not change 
the spectral properties of the graph. Any graph can be described as in section || but the price 
to pay is to add some vertices with coordinence 2, which increases V and B. In sections || and [7| 
one has already noticed that vertices of coordinence 2 play no role in the construction of orbits. 
On the other hand, since the sizes of the matrices that we have introduced are related to B and 
V, the calculations would become easier if it would be possible to minimize those numbers. It is 
the purpose of this section to explain how to generalize ( ^J36| ) to the more general case where 
a loop can be present at a vertex and several bonds link two vertices. A discussion is also given 



in [10 1 using a different method. We set 7 = 1. 

a 




Figure 9: 

(i) The expression (|34|) is again a good starting point. One considers the determinant of a graph 
one part of which is represented in figure || In (^) , the terms that involve vertex /i are: 



Since does not appear anywhere else one can integrate over <^>„, using the completeness 
relation for the propagator. One finds 

S{<i)~ x = / ■■■Gi^A^Ja)--- (155) 

In ( |155| ), the integral contains two propagators that propagate the field from vertex a to vertex 
j3: +i ((f>p,<j> a ) and Gi a0 ((f>p,<j> a ). It is now easy to see what is the natural generalization 



29 



of (35,^) if B a p bonds of lengths l 3 a g link the two vertices: 



M Ql 3 



O a 



^ coth(V7^) + 

3=1 



B, 



a/3 



e a P 



y 

fri sh(^) 



(156) 



(157) 



Let us stress that the elimination of a vertex of coordinence = 2 is particularly straight- 
forward with expression (jPD, since it is a consequence of the completeness relation for the 
propagtor of the two-dimensional harmonic oscillator. 



a 




Figure 10: 

(ii) If one now considers a graph as in figure [0], one can use the same trick as before. It is then 
easy to see that the contribution of the loop is 



L afi~''afj, v ' ' 



(158) 



where ZL, and ZjL are the lengths of the two bonds. Writing l aa = 1^ + l a for the length of 
the loop and 9 aa = 0* + 6^ a for the flux that pierces it, it is easy to see that the loop gives 
only a contribution to M aa : 



M aa = 2coth(V7U) - 2 1 C °„ ; + 



(159) 



sh(^7/ c „„, 

To conclude, gathering the points (i) and (ii) together one gives the general expression for 



M: 



M a j3 = $a(3 



y J2 coth ( + 2 yl coth(^: 



act) 



COS 



(data) 



3=1 



3=1 



sh(^/7^, 



Ee <*0 
i= i sh(^7/ 7 a/3 ) 



(160) 

where L a is the number of loops at vertex a (we did not change f\ the definition of a a p which is 
still or 1 depending the vertices are connected or not; a aa = 0). Let us remark that M a p may 



be expressed in the more condensed form jq, 10] 



M a /3 = d a p 



y coth^) - 2^ 



arcs b starting 

from vertex a 



cos(9jg) 
\ sh( % /7^ a ) 



arcs b linking 

vertices a and (3 



(161) 



where in the first sum over all arcs starting from vertex a, the case of a loop brings twice the 
same contribution, as in ( |159| ), since it is associated with two arcs starting from the vertex. 



3 this implies that a a p is no more the coordinence of vertex a. One has: m a = a a pB a (i + 2L a . 
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As an example let us consider the graph of figure ||. The number of vertices can be reduced 
to 2 and the spectral determinant is given by a determinant of a 2 x 2 matrix: 



5(7) = sh/sh6 



coth 6 + 2 coth I 
i_ 

shb 



) cos 9 
' shl 



1 

shb 
coth b 



(162) 



which gives (XT 



D Some properties of the matrices Q, R, D, e and E 

In the more general case, in the presence of a magnetic field and with mixed boundary conditions, 
the four matrices that we have introduced in section || read 

Q(a/3)(fj,v) = UapanvOaA 0p v J e 2 ; (163) 

V TTl a 1 "^y / 

R {ap){^v) = a a pa^5 av 5 PpL e~^ 1 ^ , (164) 

£{ap){ju/) = aapa^pA \7 ~ 5av) (165) 

and 

D(afl)Qi») = a a pa^5 a ^ Pv e~^ l ^ + ' ie ^ . (166) 

The inverses of those matrices are easily calculated. They are given by their hermitic conjugates 
in which one replaces y/j by —^7: Q _1 = Q*\ etc (let us recall that the transposition 

is defined as Qj a(3)(pu) = Q^u)( a (3))- F ° r example 

£ (i)(^) = a a^uSau( - <W) • ( 167 ) 

This shows that in the positive part of the spectrum, when 7 = —k 2 , the five matrices Q, R, 
D, e and E are unitary matrices. Let us notice that for Neumann boundary conditions X a = 0, 
then Q~ l = Q and e" 1 = e T . 

We now compute the determinants of those matrices. To start with, one considers the sub- 

matrices Q a defined in (f47|). We write the matrix element (Q a )f3v = ( m +A 7^/7 ~~ ^P") e ' 2 ' 
where it is understood that (3 and v belong to the set of m a neighbours of a. 

To compute the determinant of the matrix, the first step is to eliminate all the phases with 
a unitary transformation. Then one has to consider a n x n matrix of the form F a p = a — 5 a p. 
It is possible to compute the determinant of such a matrix using a recursive method, one finds 
det(F) = (-l) n (l-na). Coming back to Q a one has det(Q Q ) = "| g Tx tt ^ . It follows 

that 

If 7 = — k 2 , this is indeed a complex number of unit modulus. 

Next one would like to compute det(R). The matrix R which couples time-reversed arcs has 
2B non vanishing elements. The determinant of R is then given by the product of all elements 
of R, times the sign of the permutation V(l, 1, • • • , B, B) = (1, 1, • • • , B, B) that exchanges each 
arc with its time-reversed; it follows that 

det(R) = (-l) B e~ 2 ^ L . (169) 
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Since e 



(U* Q{~f)R{~{ = 0)U) T , it is easy to deduce its determinant: 



V 



m a - A a /^/7 
m a + A Q /^/7 



det(e) = (-l) v - B n 



(170) 



a=l 



which involves the topological invariant V — B. 

The determinant of E is given by det(£?) = det(e)e~ 2 ^ L . For Neumann boundary conditions 
det(E) = (-i) v - B e - 2 V^ L . 
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